Biosynthesis of eriodictyol in citrus waster by endowing P450BM3 activity of naringenin hydroxylation.

The flavonoid naringenin is abundantly present in pomelo peels, and the unprocessed naringenin in wastes is not friendly for the environment once discarded directly. Fortunately, the hydroxylated product of eriodictyol from naringenin exhibits remarkable antioxidant and anticancer properties. The P450s was suggested promising for the bioconversion of the flavonoids, but less naturally existed P450s show hydroxylation activity to C3' of the naringenin. By well analyzing the catalytic mechanism and the conformations of the naringenin in P450, we proposed that the intermediate Cmpd I ((porphyrin)Fe = O) is more reasonable as key conformation for the hydrolyzation, and the distance between C3'/C5' of naringenin to the O atom of CmpdI determines the hydroxylating activity for the naringenin. Thus, the "flying kite model" that gradually drags the C-H bond of the substrate to the O atom of CmpdI was put forward for rational design. With ab initio design, we successfully endowed the self-sufficient P450-BM3 hydroxylic activity to naringenin and obtained mutant M5-5, with kcat, Km, and kcat/Km values of 230.45 min-1, 310.48 µM, and 0.742 min-1 µM-1, respectively. Furthermore, the mutant M4186 was screened with kcat/Km of 4.28-fold highly improved than the reported M13. The M4186 also exhibited 62.57% yield of eriodictyol, more suitable for the industrial application. This study provided a theoretical guide for the rational design of P450s to the nonnative compounds. KEY POINTS: •The compound I is proposed as the starting point for the rational design of the P450BM3 •"Flying kite model" is proposed based on the distance between O of Cmpd I and C3'/C5' of naringenin •Mutant M15-5 with 1.6-fold of activity than M13 was obtained by ab initio modification.

A precise swaying map for how promiscuous cellobiose-2-epimerase operate bi-reaction.

Promiscuous enzymes play a crucial role in organism survival and new reaction mining. However, comprehensive mapping of the catalytic and regulatory mechanisms hasn't been well studied due to the characteristic complexity. The cellobiose 2-epimerase from Caldicellulosiruptor saccharolyticus (CsCE) with complex epimerization and isomerization was chosen to comprehensively investigate the promiscuous mechanisms. Here, the catalytic frame of ring-opening, cis-enediol mediated catalysis and ring-closing was firstly determined. To map the full view of promiscuous CE, the structure of CsCE complex with the isomerized product glucopyranosyl-ß1,4-fructose was determined. Combined with computational calculation, the promiscuity was proved a precise cooperation of the double subsites, loop rearrangement, and intermediate swaying. The flexible loop was like a gear, whose structural reshaping regulates the sway of the intermediates between the two subsites of H377-H188 and H377-H247, and thus regulates the catalytic directions. The different protonated states of cis-enediol intermediate catalyzed by H188 were the key point for the catalysis. The promiscuous enzyme tends to utilize all elements at hand to carry out the promiscuous functions.

Rational Design of P450 aMOx for Improving Anti-Markovnikov Selectivity Based on the "Butterfly" Model.

Aromatic aldehydes are important industrial raw materials mainly synthesized by anti-Markovnikov (AM) oxidation of corresponding aromatic olefins. The AM product selectivity remains a big challenge. P450 aMOx is the first reported enzyme that could catalyze AM oxidation of aromatic olefins. Here, we reported a rational design strategy based on the "butterfly" model of the active site of P450 aMOx. Constrained molecular dynamic simulations and a binding energy analysis of key residuals combined with an experimental alanine scan were applied. As a result, the mutant A275G showed high AM selectivity of >99%. The results also proved that the "butterfly" model is an effective design strategy for enzymes.

Anchor-Locker Binding Mechanism of the Coronavirus Spike Protein to Human ACE2: Insights from Computational Analysis.

COVID-19 has emerged as the most serious international pandemic in early 2020 and the lack of comprehensive knowledge in the recognition and transmission mechanisms of this virus hinders the development of suitable therapeutic strategies. The specific recognition during the binding of the spike glycoprotein (S protein) of coronavirus to the angiotensin-converting enzyme 2 (ACE2) in the host cell is widely considered the first step of infection. However, detailed insights on the underlying mechanism of dynamic recognition and binding of these two proteins remain unknown. In this work, molecular dynamics simulation and binding free energy calculation were carried out to systematically compare and analyze the receptor-binding domain (RBD) of six coronavirus' S proteins. We found that affinity and stability of the RBD from SARS-CoV-2 under the binding state with ACE2 are stronger than those of other coronaviruses. The solvent-accessible surface area (SASA) and binding free energy of different RBD subunits indicate an "anchor-locker" recognition mechanism involved in the binding of the S protein to ACE2. Loop 2 (Y473-F490) acts as an anchor for ACE2 recognition, and Loop 3 (G496-V503) locks ACE2 at the other nonanchoring end. Then, the charged or long-chain residues in the ß-sheet 1 (N450-F456) region reinforce this binding. The proposed binding mechanism was supported by umbrella sampling simulation of the dissociation process. The current computational study provides important theoretical insights for the development of new vaccines against SARS-CoV-2.

Identification of three new compounds that directly target human serine hydroxymethyltransferase 2.

Mitochondrial serine hydroxymethyltransferase 2 (SHMT2) is an important drug target in the one-carbon metabolic pathway, since its activity is critical for purine and pyrimidine biosynthesis. Additionally, it plays a prominent role during metabolic reprogramming of cancer cells, and SHMT2 inhibitors have proven useful as anticancer drugs. Compared to drugs targeting one-carbon metabolic enzymes (mainly dihydrofolate reductase and thymidylate synthase) that have been used for clinical treatment of cancer, efficient SHMT2-specific inhibitors are lacking. Therefore, we established a direct system for virtual screening, protein expression, and identification of inhibitors targeting SHMT2. First, 27 compounds qualifying as potential SHMT2 inhibitors were selected for biological activity verification through virtual screening of the 210 thousand compounds registered in the Specs database. Second, these 27 hits were subjected to quick screening by an in vitro non-competitive kinetic assay of SHMT2 single-enzyme catalysis. This allowed us to identify three compounds featuring medium-strength and non-competitive inhibition of SHMT2: AM-807/42004511 (IC50  = 14.52 ± 4.1665 µM), AM-807/40675298 (IC50  = 12.74 ± 5.8991 µM), and AM-807/42004633 (IC50  = 9.43 ± 0.5646 µM). We describe a quick screening method for the identification of inhibitors targeting SHMT2, providing a basis for subsequent identification and screening of new inhibitors.

Engineering the biomimetic cofactors of NMNH for cytochrome P450 BM3 based on binding conformation refinement.

Cytochrome P450 BM3 (BM3) is an important oxidoreductase that is widely used in drug synthesis, chemical synthesis, and other industries. However, as BM3 unquestionably increases costs by consuming a natural cofactor that unstably provides electrons, an alternative biomimetic cofactor with simpler structures represented by nicotinamide mononucleotide (NMNH) has been utilized. Currently, few reports exist on artificially modified BM3 enzymes using NMNH, especially regarding theoretical simulation and calculation. With the cognition of the mechanism in mind, we propose a strategy that optimizes and refines catalytic conformation. Based on constrained molecular dynamics simulation, the distance between N-5 of FAD flavin and C-4 of NMNH is used as a cue for the determination of improved conformation, and the potential positive mutants are subsequently screened virtually in accordance with binding free energy requirements. As a result, the K cat/K M values of the favorable mutant S848R increased to 205.38% compared to the wild-type BM3 with NMNH. These data indicate that our strategy can be applied for the specific utilization of biomimetic cofactors by oxidoreductases represented by BM3.